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INTRODUCTION 


Many current composites are made with rather brittle thermoset resins and 
have low interlaminar fracture toughness. As a result, these laminates are 
easily damaged. An understanding of the interaction between the fibers and 
the resin during interlaminar fracture could provide useful guidelines for 
developing tougher composite systems. The purpose of this study was to 
contribute to this understanding through an analysis of the delaraination front 
in a double cantilever beam (DCB) specimen. The DCB specimen has been widely 
used to characterize the mode I interlaminar fracture toughness of laminates. 
Several experimental [1-6] and analytical [7-11] studies have been conducted 
to analyze the fracture mechanism and factors influencing fracture. The 
effects of adherend configuration and material properties on the stress 
distribution and the strain-energy-release rate were investigated in 
references 5, 8, and 10. The influence of adhesive thickness on the amount of 
yielding ahead of the delamination front was also explored in reference 10. 

In the analyses reported to date, the DCB specimen has been modeled as 
two homogeneous, orthotropic adherends with or without a resin interface 
layer. No attempt has been made to examine the stress state within the 
adherend by modeling the fibers and resin separately. The present study 
addresses this need by analyzing a DCB specimen using a fiber and resin 
micromechanics model of a small region at the delamination front. A three- 
dimensional (3D) model [11] of the complete DCB specimen with homogeneous 
material properties was analyzed to determine the boundary conditions for the 
local fiber-resin model . The present study had the following objectives: (1) 

to model the delamination region of a DCB specimen, representing discrete 
fibers and resin, (2) to analyze the stresses within this fiber-resin model, 
and (3) to estimate the extent of yielding in the resin. 



The local model had a height of one ply thickness and extended a little 
more than one ply thickness ahead of and behind the delamination front. The 
fiber-resin portion of the local model contained four fibers and surrounding 
resin, with a thin “resin-rich" interface layer, typical of cocured 
graphite/epoxy laminates. The remainder of the local model was homogeneous. 

A finite element analysis with twenty-noded, 3D, parabolic elements was used. 
The displacements calculated from the 3D analysis of an orthotropic DCB 
specimen [11] were imposed as boundary conditions on the local model. 

Stress components within the fiber-resin region were calculated and 
stress distributions are presented along selected planes and surfaces within 
the local fiber-resin model. Yielding in the resin interface layer and the 
fiber-resin region was estimated using computed elastic stresses with the von 
Mises and a modified von Mises yield criteria. The modified von Mises yield 
criterion accounted for hydrostatic stress effects. The extent of yielding in 
the composite DCB specimen was also compared with that from an all-resin DCB 
specimen. 
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Poisson's ratio of resin 
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ANALYSIS 

The 3D analysis of a DC specimen [11] showed that the stress state was 
nearly uniform along the delamination front except near the specimen edges. 
As a result, the delamination stresses can be studied by examining a typical 
interior ’’slice' 1 of the specimen. Details of the analysis are presented in 
this section. 


Specimen Configuration and Materials 
Figure 1(a) shows the DCB specimen consisting of two cocured adherends 
with a resin interface layer 10 /xm thick. Each adherend represents a 12 -ply 
unidirectional graphite/epoxy laminate 1.65 mm thick. The delamination, 
located in the middle of the resin interface layer, has a length a of 50.8 
mm and width w of 25,4 mm. The elastic properties used for the fiber, 
resin, and graphite/epoxy lamina are given in Table 1. The resin and lamina 
properties were taken from reference 11 and used in micromechanics equations 
[12] to calculate the fiber properties by iteration. The specimen is loaded 
as shown, by imposing uniform displacements in the y-direction. The 
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corresponding load per unit width is denoted as P. Unless otherwise 
specified, the results in this report correspond to P =* 1 N/m. 

Local Region Modeling 

Figure 1(a) also shows a small 3D, rectangular region at the delamination 
front. This small region represents a typical slice of the DCB interior. As 
shown in figure 1(b), this local region is sub-divided into a fiber-resin 
region and a homogeneous orthotropic region. The elastic properties of the 
homogeneous region were the same as used in the global DCB model. The local 
model has a height of about one ply thickness (135 fim) and extends 165 jum 
behind and ahead of the delamination front. 

For simplicity, the fibers were arranged in a regular square array [13] 
with a fiber volume fraction of 0.63, typical of graphite/epoxy laminates. 

For a typical fiber diameter of 7 /un, the center distance between adjacent 
fibers was 8 ^m. Figure 2(a) shows the arrangement of fibers in the fiber- 
resin region. Due to symmetry, only half of each fiber was modeled, giving a 
model thickness of 4 jum. Computational limitations restricted the number of 
fibers to four. Thus, the fiber-resin region was 37 /im high and extended 80 
/im behind and 165 /im ahead of the delamination front. For convenience, the 
same finite element mesh was used for all x « constant planes (including the 
homogeneous, orthotropic portion of the local model). Also, the mesh 
refinement behind the delamination front was a mirror image of the mesh ahead 
of the front. The mesh refinement at the delamination front is shown in 
figure 2(b). In the collapsed elements at the delamination front, the mid- 
side nodes were moved to the quarter points [14, 15] to produce a stress field 
with a square root singularity. The size of the collapsed elements was 0.05 



Several analyses were performed by varying the mesh refinement in the x- 
direction and in the y-z plane. Stresses along the delamination plane were 
compared with each other and with the fine-mesh, 2D results from reference 10. 
A 3D, coarse mesh with 1892 elements and 10,670 nodes was found to adequately 
describe the singular stress distribution ahead of the delamination. However, 
this coarse mesh did not give smooth stress distributions in the width (z) 
direction. This problem was solved by refining the mesh in the y-z plane to 
obtain the model used throughout this study. This model, shown in figure 2, 
had 4578 elements and 23,037 nodes with 69,111 degrees of freedom. 

Boundary Conditions for the Local Model 
The stress analysis of the local model was performed by imposing nodal 
displacements on the faces of the model (except at the delamination plane). 

As previously mentioned, these displacements were calculated from the 3D 
analysis of the DCB specimen in reference 11. On the back face of the model 
(z * 0) , the w-displacements were set equal to zero. On the delamination 
plane (y * 0) , the stresses were zero behind the delamination front and the v- 
displacements were zero ahead of it. A brief description of the other 
boundary displacement conditions is given in Appendix A. 

Strain-Energy-Release Rate Computation 
The crack- opening- displacement method was used to calculate the stress- 
intensity factor Kj and then the strain- energy- release rate was 

calculated from K^. . Assuming the stress state at the delamination front was 
nearly plane strain, the equation for the stress - intensity factor is given by 
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where v- is the opening displacement at the first quarter-point node (8 - 
0.0125 /im) behind the delamination front. This is one -half the crack 

opening displacement. The constants and represent the resin layer 

elastic modulus and Poisson's ratio, respectively. The strain-energy-release 
rate G^. is obtained from using the following equation 


G i = 


2 2 
(1 - OK f 

r I 


( 2 ) 


RESULTS AND DISCUSSION 


The formulation of the local model was evaluated by performing an 

analysis of the local region, assuming it to be homogeneous and orthotropic, 

and then comparing these results with the 3D analysis of the DCB specimen 

[11]. The calculated stress distributions ahead of the delamination front 

from the two analyses were nearly identical . The value of G ^ calculated from 

-4 2 

the local analysis was 0.592 x 10 J/m , which agreed very well with 0.597 x 
10' 4 J/m^ for the mid- section of the delamination front from reference 11. 
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Stresses in the Interlaminar Resin Layer 

Figure 3, taken from reference 11, shows the distributions of the o^, 

0 ^, and stresses along the x-axis in the resin layer ahead of the 

delamination front. This figure shows the very high stresses and stress 
gradients immediately ahead of the delamination front. Of the three stress 
components, a is dominant; hence, the following discussions will focus on 
this stress. 

Figure 4 shows the o distribution ahead of the delamination front 
calculated from the local region analysis. A logarithmic scale was chosen to 
emphasize the region close to the delamination front. The stress 
distributions at z - 0 and z - 4 pm are virtually identical. Both stress 
distributions show the characteristic slope of -1/2 very near delamination 
front and have shapes similar to that of curves from the 2D analyses in 
references 7 and 10. Hence, the fibers modeled discretely in the present 
analysis did not influence the stresses along the x-axis ahead of the 
delamination. 

Figure 5 shows a ^ stress distributions in the resin layer immediately 
above the delamination front. These distributions are in the x * 0 plane 
along lines parallel to the delamination front. Very close to the 
delamination front, y * 0.75 pm, the a stress is nearly constant over the 
model thickness. Farther away from the delamination, the a stresses are 
lower, as expected, but are slightly elevated near the fiber centerline (z - 4 
pm ) . The load path through the fibers is stiffer than along a parallel path 
through the resin; therefore, the a stress is elevated under the fiber. 
However, as figure 5 shows, this trend dissipated within the resin layer. As 
a result, the o stress is virtually constant along the delamination front 
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over the interior portion of the DCB specimen represented by this fiber-resin 
model . 


Strain- Energy-Release Rate 

Figure 6 compares the values of the strain- energy- release rate G I along 

the delamination front for three models of the DCB specimen: (1) a 2D plane 

strain model [10], (2) a 3D homogeneous, orthotropic model [11], and (3) the 

present fiber-resin model. For the fiber-resin model, at the fiber 

centerline is only about two percent higher than at z = 0 . This 

distribution is nearly constant because a ^ varied very little along the 

delamination front, as shown in figure 5. The average for the fiber- 

-4 2 

resin model is 0.564 x 10 J/m . This is about seven percent lower than the 

-4 2 

3D value of 0.597 x 10 J/m for the specimen midplane [11]. The average 
value for the fiber-resin model is only about four percent lower than the 2D 
plane-strain value of 0.570 x 10 ^J/m^ [10]. 


Ply Stresses 

Figure 7 shows the a ^ stress versus z through the fiber-resin model. 
The solid curve, from figure 5, is the stress distribution at the "interface” 
between the ply and the resin layer. The other three curves represent the o ^ 
distributions midway between the fibers. For each curve, the a ^ stress is 
highest where the fibers are closest together, at the fiber centerline (z « 4 
pm). Also, as expected, the o ^ stresses decrease as y increases away from 
the delamination. 


Figure 8(a) shows resin stress distributions in the x - 0 plane. 

The dashed curves represent the o ^ stresses for z = 0 (the y-axis) and the 
solid curve segments represent stresses for z * 4 pm (the centerline 
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through the fibers). For comparison, the stress distribution for an all- 

resin model is also shown. For z * 0, the dashed a curve varies in a 
cyclic manner and is lower than the o ^ curve for the all -resin model. 
However, the solid curve segments for z * 4 /*m are higher than the all- resin 
(dash- dot) curve, except near the delamination front. The fibers produce 
stress concentrations in the resin where the fibers are closest together. 
Figure 8(b) shows similar a ■ stress distributions . Again, the solid curve 

segments for z * 4 /im are much higher than the dash-dot curve for the all- 
resin case, indicating a significant a ^ stress concentration between the 
fibers. Figure 8(c) shows similar stress concentrations between the 

fibers. These local stress concentrations will be discussed later in terms of 
their influence on resin yielding near the delamination front. 

Fiber-Resin Interface Stresses 

To examine the stresses at the fiber-resin interfaces, cylindrical 

coordinates were used for each fiber. As shown in figure 9, r, 9 , and x 

represent the radial, circumferential, and axial fiber directions. The 

interfacial stress state consists of a radial stress and two shear 

stresses a A and a , as shown. Figure 9 shows the interface stress 
r0 rx & 

distributions for the first fiber above the delamination front (in the x - 0° 
plane). This is the most highly stressed portion of the fiber-resin 
interface. As expected, the stress has peak values at 0 - 0° and 180° 

and a minimum at 0 - 90°. Because of symmetry, a ^ is zero at 9 - 0° and 
180°. This dashed curve changes sign near 6 = 90° and has peaks near 45° and 
160°. The distribution has peaks at 0° and 180° and changes sign near 

0-90°. 
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Next, each of the three curves discussed in figure 9 is compared with its 
counterparts for the other three fibers. Figures 10(a), (b) and (c) compare 
the , and stresses, respectively, for the four fibers. These 

curves have similar shapes, with the expected lower magnitudes for fibers 
farther from the delamination front. 

Distributions of the interface stresses a and a at 9 = 0 along 

the first fiber are shown in Figure 11. The normal stress a ^ reaches a peak 
just ahead of the delamination front (at x * 0) and then rapidly decreases. 
The shear stress has its maximum value slightly behind the delamination 

front and changes sign ahead of the delamination front before decreasing to 
zero. Since the fiber-resin interfaces near the delamination front are 
subjected to combined normal and shear stresses, interfacial strength and 
toughness analyses will, therefore, probably require multi-axial stress 
criteria. 


Yielding Near the Delamination Front 

The computed elastic stresses were used with the von Mises yield 

criterion to estimate the region of resin yielding near the delamination 

front. The stresses used in calculating the yield zone were found by scaling 

the load until G T was equal to G and, therefore, corresponded to the 
1 it 

incipient delamination growth condition. However, because this procedure 
uses elastic stresses and does not account for stress redistribution due to 
yielding, the resulting yield zone estimates should be smaller than actual 
values. The yield zone for the fiber-resin model was compared with that for a 
all-resin DCB specimen loaded to the same G^ condition. Also, yield zones 
were compared using the von Mises yield criterion and a modified von Mises 
criterion that includes hydrostatic stress effects. 
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Figure 12 shows calculated yield zones based on the von Mises yield 

2 

criterion for a load level corresponding to - 85 J/m (a typical G^ c 

a brittle graphite/epoxy composite [2]). This rather low value of G^ c was 

selected so that yielding did not develop beyond the fiber-resin portion of 

the local model. Figure 12(a) shows a cross-section at x - 0 and 12(b) shows 

the front surface of the model (z - 4 pm ) . The expected yielding in the 

interlaminar resin layer ahead of the delamination was found. However, 

localized yielding was also found between adjacent fibers. As indicated 

previously (figure 8), these regions have resin stress concentrations. A 2D, 

finite element analysis of an all-resin DCB specimen was performed for a load 

2 

corresponding to G^. - 85 J/m . The yield zone height was found to be less 

than the thickness of the interlaminar resin layer in the fiber- resin model. 

A comparison of the yield zones in the fiber-resin model and the all-resin 

model showed that the yielded volume in the fiber -res in model was about 3.5 

times that for the all-resin case. The resin stress concentrations caused by 

fibers increased yielding compared to the all-resin case. This contradicts 

the widely held assumption that fibers restrict yielding at a delamination 

and, therefore, cause a smaller yield zone than in a corresponding all-resin 

case. Of course, the present procedures provide only approximate yield zones 

2 

and the cases compared represent rather brittle materials (G^ = 85 J/m ) . 

Because yielding contributes to fracture toughness G^ , the present 
yield zone comparison may provide another explanation for the observation in 
reference 2 that G ^ for a DCB specimen can be twice as large as that for an 
all-resin specimen. This unexpected observation was explained in reference 2 
by fiber bridging effects that can elevate the interlaminar toughness compared 
to the all- resin value. The effective height of the yield zone estimated for 
the present fiber-resin case was found to be about twice that for the 
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corresponding all-resin case. Thus, the present comparison may 
be quantitatively significant; however, it is limited to rather brittle 
materials and cannot address the important trends shown in reference 2 for 
tough resins and composites. 

The yielding of thermoplastic polymers has been shown to be influenced by 
the hydrostatic component of the stress state [16-20]. In the present 
problem, the stress state near the delamination front was found to have a 
significant hydrostatic tension component. For thermoplastic resins this can 
increase yielding compared to that predicted by the von Mises criterion. 
Although the present analysis used elastic material properties and a G^ c 
that represented a thermoset resin, the computed stress distributions should 
be similar to those for a thermoplastic resin. Therefore, a modified von 
Mises criterion (Appendix B) , which includes the hydrostatic stress component, 

was used to obtain a second estimate the yield zone. This calculated yield 

2 

zone, again corresponding to « 85 J/m, is shown in figure 13. A 
comparison of figures 12 and 13 shows a significant increase in the size of 
the yield zone due to the hydrostatic stress effect. 

CONCLUDING REMARKS 

A 3D finite element model was developed to analyze the fiber -res in 
interaction at the delamination front in a double cantilever beam (DCB) 
specimen. This model represented a small portion of a graphite/epoxy ply as 
well as one-half of the thin resin interface layer with the delamination at 
the DCB specimen midplane. Within this model, the fibers were arranged in a 
square array with a fiber volume fraction of 0.63. Resin stress 
distributions, fiber-resin interface stress distributions, strain-energy- 
release rates, and resin yielding were analyzed. 

12 
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1. Within most of the interlaminar resin layer, the delamination opening 
mode stress a ^ was nearly uniform across the specimen interior. This a ^ 
stress was slightly higher in the region where the fiber was closest to the 
delamination. 

2. The fiber had a small influence on the delamination strain-energy- 

release rate , which varied by only about two percent because of the fiber. 
The average value for the fiber -res in model was about four percent 

smaller than the 2D, plain strain value for an equivalent homogeneous, 
orthotropic DCB specimen. 

3. Near the delamination front, the fiber-resin interfaces were subjected 
to combined normal and shear stresses that varied around and along the fibers. 
These results suggest that a multi-axial stress criterion may be required to 
analyze the fiber-resin interface strength or toughness. 

4. Resin stress concentrations were found between the fibers. These 

stress concentrations produced localized, fiber- induced yield zones between 
the fibers. The volume of the estimated equivalent yield zone for this fiber- 
resin model was larger than that for a corresponding all -resin DCB model 
loaded to the same level. This suggests that the fibers can increase the 

extent of matrix yielding associated with delamination growth and may, 
therefore, increase fracture toughness, rather than restrict it as usually 
assumed. 

5. The stress analysis indicated a significant hydrostatic tensile stress 
near the delamination front. When a modified von Mises yield criterion was 
used to account for the hydrostatic stress, larger yield zones were estimated 
for the delamination front. Thus the effects of hydrostatic stress on 
yielding onset and subsequent inelastic deformation should be included in 
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future DCB specimen analyses. This should be especially important for resins 
that craze under the influence of hydrostatic tensile streses. 
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APPENDIX A.- DISPLACEMENT BOUNDARY CONDITIONS FOR THE LOCAL MODEL 


This appendix describes the displacement boundary conditions used with 
the local finite element model shown in figure 1(b). This local model 
represents a typical slice of the interior of the DCB specimen in figure 1(a). 
The boundary displacements for the local model were calculated from a 3D 
finite element analysis of this DCB specimen. The fiber and resin properties 
(with a fiber volume ratio of 0.63) used in the fiber-resin model were 
equivalent to the homogeneous, orthotropic properties used in the 3D 
analysis. Details of this global analysis are given in reference 11. Only 
the displacements needed as boundary condition for the local model are 
presented here. 

The u- and v-displacements varied very little across the middle half of 
the 3D model (less than 0.05 percent for u-displacements and 0.4 percent for 
v-displacements). Therefore, these displacements were assumed to be constant 
over the local model thickness. Figures 14(a) and (b) show the distributions 
of u and v, respectively, for the top edge of the local model. Each of 

these displacement distributions was fit by a cubic spline interpolation. The 

resulting equations were then evaluated for each finite element node on the 
top edge to establish the nodal displacement. The u- and v-displacements for 

the ends of the model are shown in figures 14(c) and (d) . On the bottom edge 

of the model, v was zero ahead of the delamination front; behind the 
delamination front, the surface was stress free. Because of the symmetry of 
the square array of fibers, u- and v-displacement boundary conditions were not 
needed for the back face (z - 0) or the front face (z - 4 /im) of the local 
model . 


i 
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Because of symmetry, the w-displacement on the back face was zero. The 

w- displacements on the front face were found to vary by less than six percent 

_3 

from an average value of -0.81 x 10 pm found from the 3D DCB analysis. 
Furthermore, because this value was three orders of magnitude less than the 
u-displacements and two orders of magnitude less than the v-displacements , 
this average w-displacement was imposed on the front face. 

Because the w-displacements were nearly uniform over the local model, its 
3D stress state could be approximated as generalized plain strain. As a 
result, the local model could have been analyzed using generalized plain 
strain boundary conditions without significant error. 



APPENDIX B - A MODIFIED VON MISES YIELD CRITERION 


The yielding of materials under multi-axial stress can be predicted by 
the following form of the von Mises yield criterion. 


3 

JT r oct 


(Bl) 


ys 

Here r _ is the octahedral shear stress and a is the uniaxial yield 
oct ys 

stress in tension. Equation (Bl) was developed from a distortional energy 
analysis and assumes that the hydrostatic stress has no influence on 

yielding. For metallic materials this appears to be reasonably accurate. 
However, for thermoplastic polymers, yielding is sensitive to the hydrostatic 
stress [16-20]. A hydrostatic compression stress increases the resistance to 
yielding, while a hydrostatic tensile stress decreases the resistance. 
Sternstein and Ongchin [18] suggested a modification to equation (Bl) to 
include the effects of the hydrostatic stress. The modified equation is 
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. + a 0 a 
oct 2m 
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(B2) 


This equation is referred to herein as the modified von Mises yield criterion. 
The constants a^ and a 2 were obtained by fitting equation (B2) to the test 
data reported in references 17 through 20 for various thermoplastic resins. 
Figure 15 shows a plot of the octahedral shear stress (normalized by the 
uniaxial yield stress in tension) as a function of the normalized hydrostatic 
stress for four polymers: polyvinyl chloride (PVC) , polystyrene (PS), 

polycarbonate (PC) , and polymethyl methacrylate (PMMA) . Using a linear 
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regression analysis, the constants a^ and a^ were found to be 1.91 and 

0.30, respectively. The resulting equation is shown in figure 15 as a solid 

line. If the hydrostatic component had no effect on yielding, the line would 

be horizontal. The negative slope shows that tensile hydrostatic stresses 

tend to reduce the octahedral stress corresponding to yielding. Equation (B2) 

reduces to equation (Bl) for the uniaxial case, where a - o / 3. 

m ys 
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Table 1.- Elastic material properties for graphite/epoxy. 


Elastic Constants 

Fiber 

Resin [11] 

Lamina [11] 

E n> GPa 

211 

3.4 

134 

E 22 * 

42.0 

3.4 

13.0 

E 33 , GPa 

42.0 

3.4 

13.0 

G^ 2 » GPa 

120 

1.3 

6.4 

G ^ ^ > GPa 

120 

1.3 

6 . 4 

G 23 , GPa 

14.5 

1.3 

4.8 

"12 

0.36 

0.30 

0.34 

"13 

0.36 

0.30 

0.34 


0.45 

0.30 

0.35 

23 





a 1,2, and 3 refer to fiber, transverse, and thickness 
directions , respectively. 
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Figure 1.- DCB specimen configuration, loading, and local model (all 
dimensions in /im) . 
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Figure 3.- Stress distributions along the x-axis ahead of the delamination fo 
the 3D orthotropic, homogeneous model [11]. 
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Figure 5.- Distributions of o in the resin layer above the delamination 




3D homogeneous, orthotropic DCB [l l] 
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Figure 6.- Comparison of strain-energy-release rates along the delamination 
front . 
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Figure 7.- Distributions of o between the fibers in the z -direction . 
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Figure 8.- Ply stress distributions in the y-direction above the delamination 
front . 
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(c) a stress distributions 
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Figure 9.- Fiber-resin interface stresses in the first fiber above the 
delamination front (in the x - 0 plane). 
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stress distributions 
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Figure 11.- Fiber-resin interface stresses along the first fiber. 
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(a) x = 0 plane (b) z = 4 pm plane (front surface of model) 





( 1 - 91 t oc ,+ 0.30 a m )/ a, 
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Figure 15.- Effects of hydrostatic stress on yielding of thermoplastic resi 
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